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ABSTRACT 


The sudden onset of the northern summer monsoon has led to 
MamveeneOries as to gthe cause, one of these being a possible 
relation to the elevated Tibetan highland that dominates the 
area. 

Using a zonally symmetric numerical model driven by a 
specified equilibrium heating function, two experiments are 
Conducted. fhe first experiment integrates the model over a 
180-day period on a flat earth, linearly changing the equilib- 
rium temperature from January to July. The second experiment 
repeats the first except for the inclusion of the topography 
Bf 85°E longitude. 

In both experiments, actual velocities are underestimated 
with the topography experiment providing the better results. 
Development of the monsoon was poor in both experiments, due 
mainly to the inability of the model's vertical resolution to 
Simulate the extreme change in topography elevation and the 
unrealistic approach of a linear temperature change. Possible 


improvements for future studies are suggested. 
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i LN ERODUGIEON 


Monsoons, meaning ''season'' from the Arabic, are complex, 
seasonally dependent weather systems covering the entire 
tropics and subtropics. The northern summer monsoon begins 
in June and is characterized by a surface monsoon trough 
@ver Imdia, with an upper level anticyclone over the Tibetan 
plateau. To the south of this Tibetan anticyclone, near 10°N 
to 15°N, is a strong easterly jet stream which encompasses 
almost the entire tropics, the maximum of the jet stream 
velocities being centered at approximately 85°E longitude. 

Of interest to meteorologists for many years has been the 
rapidity with which the onset of the monsoon activity occurs. 
During the winter months, the Asia continent is under the 
influence of a subtropical westerly jet stream at the upper 
levels. Over avery brief period of time, measured in days, 
the westerly jet stream migrates toward more northern lati- 
tudes and is replaced in the tropics by the tropical easterly 
jet stream. Coupled with this replacement is the surface 
monsoon trough development and its associated surface south- 
westerly winds, winds laden with moisture from the oceans. 
These southwesterly winds cause widespread rainfall as they 
reach the southern Asia continent. 

The need to understand the sudden onset of the summer mon- 
soon has led to many theories among meteorologists. Yin (1949) 


has proposed that this onset, as represented by the sudden 








BemmaccMentOL the Wwesaterty jet by the. easterly jet in the 
BEepica, May be related in some way to the effect of the 
iopetcan haghlands, an elevation that is very striking in 
Mac Ure. 

The purpose of this study is to examine such a possibility 
by using a numerical model to simulate the development of the 
summer monsoon with different experiments. For the first 
experiment, a zonally symmetric January steady state condition 
mime tliat Garth 1S obtained. Starting from this condition, 
the model is integrated for 180 days with the mean thermal 
condition changing from January to July. The final state 
represents a July steady state condition. In the second ex- 
periment, the distribution of topography along 85°E, which is 
fear the center longitude of the Tibetan highland, is included 
Paeene model. The terdagy integration is repeated. 

It is hoped that by comparing the results of these two 
experiments, some understanding of topographical effects upon 
the evolution of events that leads to the development of the 


Summer monsoon may be obtained. 





IT. MQDEL AND PROCEDURE 


fee 6d BAAN LC CMODEL 

In this study, the model formulated by Monaco and Williams 
icy 5) as used. This Wedel uses the primitive squeticne in 
ppictical Coordinates with ‘sigma” as the vertical coordinate. 
ine homazental distribution of variables is on a 4-degree 
resolution. The model is made zonally symmetric by removing 
mebezonal Varilattons. The troposphere 1s divided into seven 
equally spaced "sigma" levels with three main reporting levels. 
The sigma coordinate is defined as: 

P ~ Py 


2 (2.1) 
TT 


where p is the pressure of the sigma level, P, 1s the height 
of the constant tropopause which is set at 100 mb, and 7 is 


mie terrain pressure defined as: 


i= ePe Pe (2.2) 


where pe 1s the surface pressure. 

Time integration is comprised of continuous sections of 
one Matsuno and four leapfrog steps with time increments of 
eemnutes. All data for this study are interpolated to con- 


meme Pressure sugraces of 250 mb, 500 mb, and 850 mb. 


10 





B. SPECIFICATION OF HEATING AND TOPOGRAPHY 

The model is driven by a heating function, H, which con- 
Gimucusly adjusts the model atmosphere temperature to a 
specified equilibrium temperature, T*, using the following 


formula: 


f= - eee te) |, iS 


where T is the model atmosphere temperature and t is a speci- 
fied adjustment time, set at two days for this study. 

Pim TOneThaoay —experimen:, the equilibrium tempera- 
tures for January and July are extracted from observations of 
Palmen and Newton [1969] for the January and July climatology 
(Figures 1 and 2, respectively). The dashed lines indicate 
the horizontal sigma levels used in obtaining the tempera- 
tures, in this case assuming the model levels to be 250 mb, 
soo Mb, and 850 mb. 

Because the temperatures in July at 85°E are considerably 
fener than the zonally averaged temperatures, the equilibrium 
temperature for July is modified by the addition of a departure 
from the zonal mean, at 85°E, as observed by Krishnamurti [1971] 
and shown in Figure 3. At the 250 mb level, the full depar- 
mee 1S added to the mean temperature. At the 550 mb level, 
Gimechcent ©: the departure is added. At the 850 mb level, 
mepercent Of the departure 1s added to the mean temperature. 

In the topography experiment, topography as a function of 
latitude is shown in Figure 4. It is introduced through a 
gradual change in the surface geopotential, ¢_, as given by 


S 
the following formula: 


Joi 
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mieure 4: Topography elevation at 85°E longitude. 


eS 


(m 


HEIGHT 





Tt 
5 7% (sin® (y76)) (2.4) 


5 


Here, $d¢ 1s the specified final surface geopotential value; 
ty is the total amount of time required for mountain building, 
specified as 24 hrs for this study; and t, is the number of 
incremental changes to occur within the mountain building 
item, specified as 48 hrs for this experiment. 

The gradual build-up of the topography is to avoid the 
undesirable numerical shock of a sudden change invariable 
magnitudes. The same January equilibrium temperature for the 
no-topography experiment is used in this experiment before 
[meen introduction Of topography, but after the topography is 
fully developed the temperature specified along the dotted 
lines in Figure 1 is used. These curves, which are initial 
Sigma levels, are patterned after the topography. 

The July equilibrium temperature for this experiment is 
extracted from Figure 2, following the same dotted lines. 
Modification of the July temperature is again accomplished 
by adding the same deviation from the mean as used in the 


no-topography experiment. 


Ge. NO-TOPOGRAPHY EXPERIMENT 

In this experiment, pressure levels are initialized at 
Standard atmospheric values with all velocities being zero. 
The initial temperature is specified as the January equilibrium 
temperature as shown in Figure 1. The model is started and 
aerowed tO run for a period of 90 days, sufficiently long to 


establish a quasi-steady state condition. 


16 





The integration is continued from this steady state con- 
dition and allowed to run for a period of another 180 days. 
mireouehnout this run, the equilibrium temperature is continu- 
Sucivechanged, by 2 Simple linear interpolation, to reach the 
equilibrium temperature of July. 

At the end of the 180-day run, the integration is allowed 
to continue for an additional 90 days with the July equilibrium 
temperature remaining fixed in order to achieve a quasi-steady 
Sate condition for July. 

There is no significant difference between the circulation 
Daeeerns at the end of the 180-day run and those at the end of 
the additional 90 days using constant July temperatures. Thus, 
it 1S apparent that the model adjusts well to the equilibrium 
heating function and any day within a particular simulation 
period may be considered to represent an instantaneous quasi- 


steady state condition. 


eee LOPOGRAPHY EXPERIMENT 

The procedure for this part of the experiment is similar 
mo the no-topography case except for the introduction of to- 
pography. To compensate for the mass changes due to the volume 
of atmosphere change, the surface pressure is specified as 
860 mb, a figure experimentally determined to produce an over 
ocean surface pressure of 1015 mb. 

icesane wanleany equilibrium temperature used for the no- 
mepeetaphy experiment 1s held fixed and the model is integrated 
for a period of 45 days, allowing the model to adjust to the 


meetrie¢ialiy low initial surface pressure. At this point, 


Wy 





Eee cOpegraphy 1S introduced as previously discussed and the 
“wibibrium temperature 1s readjusted to conform with the 
mopoesraphy sigma levels as shown by the dotted line in Figure 
1. An additional 45 days is allowed after the introduction 
Of topography and January temperature change to achieve a 
January steady state. 

PaOimenicusecady Condition, the integration is allowed to 
proceed for 180 days, again allowing the equilibrium tempera- 
meme to Change linearly from January to July, the July temper- 
ature being obtained as previously discussed. 

The July equilibrium temperature is again held fixed and 
the model is allowed to continue for an additional 90 days. 
There is again no significant difference between the wind 
patterns at the end of the 180-day integration and the 


additional 90-day run. 


18 








Po ae keou LTS 


imprinewel Cimeulation tLeactures Of interest in this 
malay are the Libetan high and the tropical easterly jet 
stream. The simulation of the development of these features 
by the model are compared with observations at 85°E, mainly 


w@ose Of Sadler [1975]. 


wee SLMULATED CIRCULATIONS IN JANUARY AND JULY 
imo; lopography Experiment 

The January steady state condition without the inclu- 
Sion of topography is represented by the zonal wind as a 
mmetioOn of latitude and elevation in Figure 5. The zonal 
mea at 250 mb as a function of latitude is compared with 
Meee ObServed by Sadler at 200 mb in Figure 6. The observed 
meeterly flow occupies the region north of approximately 14°N 
feeenea well detined jet core at 32°N. The maximum velocity 
mp this core 18, on the average, 44 m sect, The model pro- 
duces a broad and weak westerly flow to northward of 4°N. 
There is a definite lack of a well defined, sharp jet maximun, 


1 Ze ASIN, fee ei iba 


the maximum velocity being 21.6 m sec. 
iPeerEuaes, the velocities are grossly underdeveloped. 

HicwewlyesStcadye State condition 15 represented by the 
zonal wind as a function of latitude and elevation in Figure 
fee ene zonal wind at 250 mb as a function of latitude is 


compared with Sadler's observation at 200 mb in Figure 8. 


ime replacement of the westerly jet stream by the tropical 
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Saaechismeeteis durte evident in both data. The observed 
easterly jet axis is centered at 12°N with a maximum velocity 
oa 22 Mm sei a In this experiment, the model produced a more 
mealistic result than the January state. The simulated 
meooreal casterily j6t is clearly developed with a core axis 
aeeo N. The maximum velocity in the core is 23 m sec me Tena 
close agreement to the observed value. 
Z. lopography Experiment 

The January steady state condition with the inclusion 
Baeeene topography is displayed in Figure 9 which gives the 
meal wind aS a function of latitude and elevation. The zonal 
peemee at 250 mb 1S again compared with Sadler's observation at 
200 mb in Figure 10. The underestimate of winds, a problem 
Mame ne no-topography case, is still evident but not to such 
an extreme degree. For example, at 32°N, the no-topography 
walue of 16 m sect Mseincreased to 26 m sect. eal ese iJ 
below the 44 m sec observation, however. The simulated 
pattern changes markedly from the no-topography case. The 
Simulation shows a clearly bimodal distribution with jet stream 


t at 20°N and 26 m sec - 


Sones of 25 m sec eo en the wlat ter 
coincides with the observed latitude of the westerly jet core. 
The latitude band between these two jet cores contains the 
fmeeetcal face of the model topography. 

iiemeclo~eoteadywestate Condition is shown in Figure 11 
by the zonal wind as a function of latitude and elevation. 


Figure 1Z gives the comparison of the simulated zonal wind 


at 250 mb and the observed zonal wind at 200 mb by Sadler as 
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Pmerunectron OL latitwde. In this case, the near match in the 
welocities in the latitudes of the tropical easterly jet has 
diminished somewhat. The jet core is shifted north to 20°N 
and the velocity is 2/7 m sae oe an overestimate from obser- 
vation by approximately 7 m sect. 

There are no good observational data at 85°E longitude 
with which to compare the 500 mb and 850 mb flow simulations.’ 


The patterns are consistent with the 250 mb patterns, however, 


as to placement of features and intensities of circulations. 


B. DEVELOPMENT OF THE MONSOON FLOW 

The zonal wind at 250 mb produced in the no-topography 
experiment is shown as a function of latitude and time in 
Figure 13. The zonal wind at 250 mb produced in the topography 
experiment is shown as a function of latitude and time in 
Figure 14. 

In both the topography and no-topography experiments, the 
seasonal reversal in the zonal wind at the 250 mb level is 
completed gradually over the six-month period. 

Mies tropical Gasterly jet first appears at approximately 
12°N in mid-March for the no-topography experiment. In the 
topography experiment, it is developed at approximately 20°N 
in early March, thus a shift of the jet axis to the north by 
Seeeenees Of latitude with the introduction of topography is 
produced. This displacement in the core is not observed and 
thus must be considered incorrect. In the real atmosphere, 
the appearance of the tropical easterly jet actually begins 


in mid-May and is centered at approximately 12°N. 
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PeetucsesOrethnercarly arrival of the easterly jet in both 
experiments, the simulated timing of the monsoon development 
Meuvery poor, with the no-topography case being, paradoxically, 
the better of the two, however. 

The gradual acceleration of the zonal winds for both cases 
memapparenciy related to the specified linear change of the 
equilibrium temperatures. In the real atmosphere, the temper- 
ature change, i.e., the reversal of north-south temperature 
gradient, occurs rather suddenly in mid-May, so that the 
Simulation of the change in temperature 1s unrealistic and 


1t causes the gradual reversal of the zonal flow. 
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IV. DISCUSSION 


The results of the two experiments indicate that the 
Simplatications used in the numerical model are unrealistic. 
These simplifications can be broken into two general areas, 
one related to the model itself and a second related to the 


data used. 


Pee mODEL SIMPLIFICATIONS 

For both experiments, a zonal symmetry assumption was 
imposed. It appears that an exclusion of zonally asymmetric 
eddy motions may fea meomtlmulacccnean le (IG1reulatieon pattern, 
due to the importance of nonlinear processes, as shown by 
Krishnamurti [1971] and others. A more severe problem of the 
ieee 1s perhaps the vertical resolution, as evident by the 
extremely poor simulation of the time evolution of the 250 mb 
zonal flow in the topography experiment. Only three levels 
mi the troposphere are apparently insufficient because of the 
steepness of the Himalaya mountains. As a result, the inter- 
polation of the simulated data from sigma levels to constant 


Mmeeesute Surfaces sutfers from the crude vertical resolution. 


B. DATA SIMPLIFICATIONS 

MacaGly shown in the model results, the zonal velocity is 
underestimated in January for both the topography and no- 
mepootaphy experiment, although the July velocities are in 


somewhat better agreement with observation. As stated 
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Pecviously, eve July equilibrium temperature specified in the 
model has been modified. from the mean by observations of 
Krishnamurti [1971]. No similar modification was made in 
January as only the zonal mean temperatures are used. Because 
the wind is quite geostrophic and is coupled with temperature 
Mmeaadtent through the thermal wind equation, an error in tem- 
perature will precipitate an error in winds. It is therefore 
necessary to improve the specification of the thermal field 
for January. 

The specification of seasonal temperature changes by the 
linear interpolation between January and July produces the 
obvious and unrealistic gradual and almost linear wind changes. 
It appears necessary that the time series of the actual temper- 
ature data must be used to adjust the thermal state of the 
model atmosphere. 

Only when this realistic temperature change and an adequate 
vertical resolution are used in the model can the mechanical 
effect of the mountains on the zonally symmetric circulation 


be adequately studied. 
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